Dirac gap-induced graphene quantum dot in an electrostatic potential 
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A spatially modulated Dirac gap in a graphene sheet leads to charge confinement, thus enabling a 
graphene quantum dot to be formed without the application of external electric and magnetic fields 
[Appl. Phys. Lett. 97, 243106 (2010)]. This can be achieved provided the Dirac gap has a local 
minimum in which the states become localised. In this work, the physics of such a gap-induced dot 
I is investigated in the continuum limit by solving the Dirac equation. It is shown that gap-induced 

. confined states couple to the states introduced by an electrostatic quantum well potential. Hence 

the region in which the resulting hybridized states are localised can be tuned with the potential 
strength, an effect which involves Klein tunneling. The proposed quantum dot may be used to 
probe quasi-relativistic effects in graphene, while the induced confined states may be useful for 
graphene-based nanostructures. 
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I. INTRODUCTION 



Quantum dots have been studied extensively over the last years because they allow fundamental laws of quantum 
' physics to be probed and they might also find applications in nanoelectronics as detectors and transistors. Although 
^ . a dot in a semiconductor heterostructure can be routinely formed, in monolayer graphene the situation is different. 

An ideal graphene sheet is gapless, its energy dispersion is linear, E—±vf\p\, and its charge carriers are massless Dirac 
particles. Therefore, they exhibit Klein tunneling, which allows massless particles to tunnel through any electrostatic 
potential barrier. This property excludes the possibility of fabricating a graphene dot simply by using electrostatic 
gates, as in common semiconductors. 

In particular, the states of an electrostatic graphene dot are deconfined, i.e., they have an oscillatory tail outside 
^ I the dot region, which is a direct consequence of the Klein tunneling. This essentially means that the electrons 
Q , can only spend a finite time interval inside the dot, and it has been shown theoretically that this time interval is 
O ' sensitive to the details of the potential profile and the energy of the quantum state;^ A uniform magnetic field leads 
to confined states (bound states) which decay exponentially outside the dot, provided the electrostatic potential rises 
' slowly compared to the magnetic vector potential*S,i2 Thus a graphene quantum dot can be formed with a uniform 
^ ' magnetic field. This property has also potential applications in graphene waveguides J^*^ 

OO ' It is experimentally possible to induce an energy gap in graphene's band structure, referred to as the Dirac gapi^ 
OO In this case the energy dispersion becomes E=±y/ (wfp)^ + A^, with 2 A the value of the Dirac gap. The importance 
\l , of the gap stems from the fact that it gives rise to carriers with mass (~A/wp) and thus a quantum well potential 
' induces confinement as in common semiconductor dotsjii^ This happens because the Klein tunneling is suppressed 
, for a particle with mass, as long as its energy lies in the gap. The quantum states with energies within the gap decay 
■ exponentially in the barrier region, thus trapping the electron in the dot for a theoretically infinite time. 
] ' This work considers the case of a spatially modulated Dirac gap that enables a graphene quantum dot to be 
I ■ formed without the application of external electric and magnetic fields.J^ This is feasible provided the gap has a local 
J> I minimum in which the electron states become localized. Gap-induced confined states couple with states induced by 
an electrostatic quantum well potential; then this property allows the region in which the resulting hybridized states 
are localised to be tuned with the potential strength. This type of dot can be used to probe quasi-relativistic effects in 
, graphene, related to Klein tunneling, while the confined states may be useful for applications. In order to manipulate 
the spin states of an electron confined in a quantum dot, magnetic fields are usually needed. Thus, forming a graphene 
dot without applying additional magnetic fields to control the charge states can be advantageous. This is one of the 
motivations of our proposal. Moreover, the density of states in a gap-induced dot can be made low enough since only 
states with energies in the gap are of importance. A magnetic dot does not have this property and it is more tricky 
to induce a low density of states. 

The value of the Dirac gap depends on the specific experimental technique and can range from a few to hundreds of 
meV. For example, a Dirac gap has been measured in graphene grown epitaxially on a SiC substrate>ii A gap opening 
has been demonstrated by controlling the structure of the interface between graphene and rutheniumii^ Further, it has 
been reported experimentally a spatially modulated and (buffer layer) thickness-dependent Dirac gapji^ Moreover, it 
has been suggested theoretically that local strain and/or chemical methods might also be employed to open- up and 
tune the Dirac gapj ^'^'^"^" — 
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This paper is organized as follows. In Sec.|TT]tlie quantum dot model is presented and the necessary conditions for 
confinement are derived in the presence of a magnetic vector and scalar potentials, as well as a spatially modulated 
Dirac gap. In Sec. IIIII a model dot is examined in the regime where the Dirac gap has the same spatial profile as the 
electrostatic potential. In Sec. IIVI some semi-analytical results are derived for the case of a piecewise-constant Dirac 
gap profile and the general properties of the gap-induced dot are presented. The interplay of the gap-induced dot 
with an electrostatic quantum well potential is studied in Sec. |Vl and Sec. IVII presents the basic conclusions of the 
paper. 



II. PHYSICAL MODEL 



A. Graphene quantum dot Hamiltonian 



For energies near the Dirac points (<± 1 eV) and in the continuum limit, a graphene quantum dot can be described 
by the effective 2x2 Hamiltonian 

77 = UFcr • (p + eA) + VI-^tAct^. (1) 

The Fermi velocity vp=j/h is assumed to be position-independent, where 7=646 meV nm is a band structure pa- 
rameter. Also, (T—{ax,cry), are the 2x2 Pauli operators acting on the two carbon sublattices, X is the 2x2 unit 
matrix, p——ih^——ih(dx^ dy) is the two-dimensional momentum operator, A is the vector potential that generates 
the magnetic field B=VxA, and V is the quantum well potential formed electrostatically, for example, due to a 
gate voltage. The last term in Eq. ([T]), referred to as the mass term, gives rise to an energy gap 2 A in the spectrum 
of graphene, where t=+I (r=— 1) corresponds to the K {K') valley. The Hamiltonian ([T]) assumes no intervalley 
coupling, which is a good approximation for most graphene samples. '^'^ 

For the dot model we consider that both V and A have cylindrical symmetry and the applied magnetic field is 
perpendicular to the graphene sheet, i.e., B=_Bz. Hence only the azimuthal component Ag is nonzero, therefore 

A = (0,Ae,0), Ae^^, (2) 

which generates the field B—Bor''~^; here, s—1 corresponds to a uniform magnetic field. Hamiltonian Eq. ([1]) is 
written in cylindrical coordinates {x^rcos6, y=r sinO) with the substitution 

- idy = e-'" (^dr - '-de^ , (3) 

and 

cr • A = — {ax sin — (jy cos 6)Ag . (4) 
A two-component solution to the Dirac equation 

i?* = E^, (5) 
with E being the eigenenergy, is written in the general form 

*/2Wexp(*m0) )■ ^'^> 

with r7i=0, ±1, ... being the orbital angular momentum quantum number. The radial components /i and /2 express 
amplitude probabilities on the two carbon sublattices of graphene, and they satisfy the two coupled differential 
equations 

{V + tA)/i +(^U + 7^:) /2 = Eh, (7a) 

{V-TA)f2 = Ef2. (7b) 
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The term U is due to the angular momentum and the apphed magnetic field, and is given by 

7(2m ~ 1) jeAg 
2-r + -?r- 

Applying the time-reversal symmetry operator iCay, with C the operator of complex conjugation, to Eqs. (j7ap 
and (j7bl) . it can be shown that the eigenenergies^^- satisfy the condition E{m, B,t)=E{1 — m,—B,—T). Also, 
E{V,t)~—E{—V,—t) due to electron-hole symmetry. 

B. Confinement in a Dirac gap-induced dot 

The quantum states of a Schrodinger dot, formed by an electrostatic quantum well, can be classified according to 
their energy. A state is confined when its energy relative to the bottom of the well is smaller than the well depth. 
This state has an exponential tail at a large distance from the quantum well. In the opposite regime, when the energy 
is larger than the well depth, the state is deconfined and it has an oscillatory tail. 

This behaviour is no longer valid for a Dirac dot formed in graphene, because the energy spectrum of graphene is 
unbound and Klein tunneling takes place. In particular, a massless Dirac electron can tunnel through any electrostatic 
potential barrier, as has been confirmed experimentally via transport measurements in p-n junctions. For this reason 
an electrostatic graphene dot cannot confine electrons. 

Nevertheless, confinement can be achieved in the presence of a magnetic field and/or an energy gap in graphene's 
spectrum, though for this to happen some specific conditions have to be satisfied, which are analysed in this section. 
As a general rule, confined states should decay asymptotically, i.e., at a large (r— >-oo) radial distance r, independent 
of their energy, whereas deconfined states should have an oscillatory tail. The former states are sometimes called 
bound and the latter quasi-bound. 

Below we employ a rigorous treatment to examine the dot states.^ We derive a single second-order differential 
equation for each radial component, by decoupling Eqs. ([7a|) and (l7b| . Then we identify the form of the state in the 
asymptotic region, i.e., at a large radial distance from the origin of the dot. The /2 component satisfies 

'^ + -ir)f+bir)h^O, (9) 



with the coefficients 



and 



a{r) 



Vi 



V+-E' 



U [V- -E){V+-E) U' 

with V±—V ± tA and the prime denotes differentiation with respect to r. The first derivative term in Eq. ([9]) com- 
plicates the analysis of the quantum states. For this reason we eliminate it by writing /2 in the form f2{r)=g{r)u2{r) 
and from Eq. ([9]) we derive that 

+ (^9 + «5) ^ + (9 + ag' + bg)u2 = 0. (10) 
If we choose g=exp(— J a/2dr), the first derivative term cancels and U2 satisfies the Schrodinger- like equation 

(fu2 

dp 

with the r-dependent coefficient 



+ fc^(r)u2-0, (11) 



kl^r)=b^---^. (12) 



Here g is not an oscillatory function of r, therefore /2 has the same confined or deconfined character as U2- Specifically 
Eq. pT|) suggests that U2 is confined only if is asymptotically (r— )-oo) negative. Otherwise U2 is deconfined. The 
same arguments are valid for the /i componenti^ 
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To proceed, we assume that for a large radial distance r both V and A are constant or have a power-law dependence; 
therefore in this case we derive asymptotically that 

Equation shows that the sign of A;| is tunable with Ag, V, A and energy E, and thus the same is valid for the 
character of the dot states. Further, the character of the states is independent of the choice of valley r=±l, although 
for a fixed m the energies are not the same for both valleys. 

The form of (r) in Eq. ((T3]) shows that the vector potential acts equivalently to the mass term, namely, both have 
the tendency to confine the states. This can be understood as follows. It can be seen from Eqs. ([7a| and (|7b|) that 
the radial Hamiltonian which acts on /i and /2 is 



[/ + 7i V-tA 

The similarity transformation P'^ Hj.P—T-Lr with the operator 

results in the transformed Hamiltonian 

«--l^rA-7|: V-U )■ 

Among the terms which dominate asymptotically, U and rA are interchanged in Hr with respect to H,.. However, 
asymptotically U'^jeAe/h, therefore in this regime rA and Ae should act equivalently. 

The necessary condition for the occurrence of confinement, i.e., A:| < asymptotically, means that at least one 
of Ag, A has to be nonzero, otherwise A:| > leading to deconfined states. Most importantly, Eq. ([T3|) shows that 
confined states can be induced even for Ag=Q and ^=0 everywhere, provided that E'^ — A^<0 asymptotically. As 
shown in Ref. |T2 this inequality can be satisfied when A is spatially-dependent; for example, when A is zero within 
a disc area, and nonzero outside that area, A=So. Then, as shown below, discrete energy levels in the range |i?|<(5o 
correspond to confined states with a large amplitude within the disc area. 



III. A MODEL QUANTUM DOT SYSTEM 



To demonstrate the arguments presented in the previous section, consider first a model quantum dot for which 
the electrostatic potential V, and the mass term A have a power-law spatial dependence. Equation (jl3p shows 
that if Ag=0, and V, A are unequal, then confined states occur only if {V — A) <0, so that fc2<0. In this case 
confinement is energy-independent, i.e., all the states independent of their energy are confined. In the special limit of 
V—A only energies which satisfy —EV<0 correspond to confined states, which implies that if V increases (decreases) 
asymptotically then confinement occurs only for positive (negative) energies. 

The complete physical behaviour can be demonstrated when V and A have the same r-dependence, for example 
paraboHc, V=Vor'^ and A=Aor^ with Vq, Aq > 0. Then from Eq. p3)) the resultant states are confined when Ao>Vb 
so that A:|<0, and deconfined when Ao<Vo and hence /c|>0. For Vq—Aq only positive energies correspond to confined 
states. 

To quantify this behaviour the two coupled Eqs. (j7al) and (j7b[) are discretised on a uniform grid and the resulting 
matrix eigenvalue problem is solved numerically, using a similar technique as in Ref. ^. The boundary conditions for 
the components /i and /2 which lead to a Hermitian eigenvalue problem [Eqs. (|7ap and (j7bl) ] are derived in Ref. H, 
and it can be shown that these are not modified by the presence of the mass term. So one component has to vanish 
at the origin and the second has to vanish asymptotically (at the boundary of the computational box). The same 
procedure is followed in all sections. 

Figure [T] shows quantum states for Aq constant, Ao=l /zeV nm~^, and three different choices of Vq. The states 
plotted are those which correspond to the lowest positive energy when Vb=0. When Vo=0 (<Ao) the states are 
confined and they have an exponential tail. This demonstrates that a graphene dot can be formed thanks to a 
spatially modulated Dirac gap without applying external electric or magnetic fields. When Vb=0.9Ao (<Ao) the 
states are still confined though the region within which they are localised is slightly modified by the presence of the 
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FIG. 1: (a) Quantum states for m—2, r=+l and different potential strengths Vo. The mass term is modeUed by A=Aor^, with 
Ao=l peV nm~^, and the electrostatic potential is modelled by V=Vor^. The potential strength Vo is from top to bottom: 
Vb/Ao=0, 0.9, 1.5. The vertical axes of the insets to the top and middle panels are on logarithmic scale, (b) As in (a) but for 
r=-l. 

potential. On the other hand, when Vo=1.5Ao (>Ao) the states have undergone a transition and they are deconfined 
asymptotically with an oscillatory tail. This behaviour of the states with respect to the potential strength is valid 
for all values of m and r=±l; though, as shown in Fig. [l] the relative amplitude of the two components depends on 
the choice of r. In particular, for r=+l, /i>/2 (/i</2) when E>0 {E<0), whereas the opposite is valid for r=— 1, 
/i</2 (/i>/2) when E>0 iE<0)M 

A similar confinement-deconfinement transition in the character of the dot states can also be induced for A=0, 
provided that the vector potential is nonzero (Ag^O) and has the same power-law dependence as the electrostatic 
potential^ This can happen because the first and third terms in Eq. (|13p have the same sign. This observation might 
be the key to fabricate a graphene dot with the help of a uniform magnetic field and an electrostatic potential, and 
some preliminary calculations suggest that this should be feasible in realistic dot designs formed in gated graphenei^i^i 

Figure [2] shows the energy level diagram in a range of confined states for the two valleys t—±1. The spectrum 
consists of two ladders (sets) of discrete levels indicating the existence of confined states. The energy ladders are 
separated by a gap which is large when the angular momentum is large. As Vq increases the confinement for negative 
energies becomes weak and for this reason the splitting between the discrete levels cannot be clearly resolved. The 
pattern of the two ladders is characteristic of confined states in graphene quantum dots and, for instance, can also 
occur when the mass term is replaced by a vector potential of the same spatial dependence<^ 

IV. INVESTIGATION OF A DIRAC GAP-INDUCED DOT 

In this section the basic properties of a Dirac gap-induced dot are investigated numerically. We assume that there 
is neither an electrostatic potential (V—O) nor a magnetic field (Ag—O) and thus confinement is due solely to the 
spatial modulation of the gap. Then confinement depends on energy, as can be seen from Eq. (|13|) . and it is achieved 
when E^ — A^<0, where A is the value of the mass term asymptotically. 

A. Graphene dot formed by a piecewise-constant Dirac gap 

First a regime for which analytical results can be obtained is examined. Consider the special limit where the 
spatially-dependent mass term A changes discontinuously. In this case A can be modelled by the expression 



A(r) = So&{r - R), 



(17) 
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FIG. 2: Energy levels of confined states for the two valleys: r=+l (solid lines) and t—~1 (dashed lines). The angular 
momentum number is m—2, the electrostatic potential is modelled by V=Vor'^ and the mass term is modelled by A=Aor^. 
Here Ao=l fieV nm~^, whereas Vb is tuned. 




r (nm) 

FIG. 3; Quantum states for a piecewise-constant mass term A given by Eq. (|17|) . The parameters are m=2, 7?=100 nm, and 
5o=80 meV. Results for both valleys, t=±1, are shown. The states correspond to the lowest positive eigenenergy. 



with Q the Heaviside step-function. For simplicity we assume that A(r)=0 for r<R, although in the following analysis 
a nonzero value can be introduced straightforwardly. As shown in the previous section, the two radial components /i 
and /2 satisfy a second-order differential equation, which for V—0 and Ae=0 has the form [(rA)^=A^] 

dVi 4(m-l)^-l E^-A^^ 

~n r"2 ■'1 ^ 2 — ~ ^^^^> 

(XT i^Y"^ 

- -^^f^ + — = 0. (18b) 
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If we assume a solution of the form ji=\frTi, then Ti have to satisfy the foUowing Bessel's differential equations 

-7^ H -, 5 -7^1 H ^ -7^1 = 0, (19a) 

dr'^ r ar 7^ 

J"2 + 5 ^"2 = 0. 19b 
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Setting k~E/j and (j'=-\/|i?^ — 5ol/7' ^he solutions for confined states can be written in the general form 

r<R, (20) 



and 



■^M = I aJm-i{kr) 
T2 I V aJ,n(kr) 



with c=a/|£^^ - So\/iE + t5o) and Of, /3 are constants that can be determined from the normalisation condition and 
the requirement that both components are continuous at r=R. Here Jm is an ordinary Bessel function of the first 
kind, and Km is a modified Bessel function of the second kind. These have been chosen since Jm is regular at the 
origin (^=0), while Km decays exponentially at large radial distances as needed for confined statesi^ The coefficients 
(a, /3, c) in Eqs. (|20|) and ([2T|) are introduced in order /i and /2 to satisfy the coupled Eqs. ([7a| and (fTb]) . 
Both components have to be continuous at r—R, leading to the following algebraic equation 

cJm-i{kR)Km{qR) = Jm{kR)Km-i{qR). (22) 

This is solved numerically, with bisection, to give the energies of the confined states in the range |i?|<(5o- Thanks to the 
properties of the Bessel functions: J -m = {—^)™ Jm and K^m—Km, the energies satisfy the condition E{m,T)=E{l ~ 
m, — r) as time-reversal symmetry requires when Ag—0. Figure [3] shows one example of the two confined components 
for the lowest positive energy. When this is much smaller than Sq the states decay quickly in the region r>R, while 
the states exhibit a kink at r—R because their first derivative are discontinuous at that point. 

B. Graphene dot formed by a continuous Dirac gap profile 

In any realistic graphene system the spatial modulation of the Dirac gap will not be perfectly sharp. For this reason 
we make a more realistic choice of the mass term, than that in Eq. (|17|) . and use the expression 

r 0, r<R 

Air) = { ^ (23) 

[-(5ocosh(^) " + <5o, R <r, 

so that asymptotically Awf^o- For all the calculations we choose for the parameters i?=250 nm and d—15Q nm, 
although these choices do not affect the main conclusions. In the limit d—>-0 we recover Eq. ([TTl) . 

Figure 2] shows the energy levels as a function of the asymptotic value of the mass term , giving rise to an energy 
gap of 2So- The (quasi) continuum of levels, corresponding to the black area in Figs.Hl^a) and (b), refiect deconfined 
states. The finite system size needed for the numerical calculations generates a gap at do=0, which increases with |m|. 
However, the system size is chosen large enough and hence it does not affect the physics of the discrete levels of the 
dot that we are interested in. These levels emerging through the continuum reflect the formation of confined states. 
The transition from continuum to discrete levels, and vice-versa, is accompanied by the appearance of anticrossing 
points. As seen in Figs. |3Ua) and (b), the discrete levels form two ladders (sets) of energy with either positive or 
negative values, separated by a gap pertinent to the angular motion. The typical splitting of the levels increases with 
So because the confinement becomes stronger. Further, the two energy ladders arc formed after a critical value of Sq 
depending on the angular momentum. This happens because the angular motion tends to delocalize the states near 
the origin of the dot, see Eq. and therefore for a given Sq not all m values produce confined states. For this reason 
there are no discrete levels for (5o<3 meV and 28 meV in Fig. |4l when m=2 and m=15, respectively. 

If \m\ is large, then i5o has to be large for the formation of confined states. This effect has a direct signature in 
the density of states (DOS) which is shown in the contour plot of Fig. SJc). For simplicity only t=+1 is considered, 
since inclusion of r=— 1 simply doubles the number of states. Moreover, because V=A0=O the energies satisfy the 
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FIG. 4: (Color online) (a), (b) Energy levels, for t—+1, as a function of the asymptotic value of the mass term So which 
generates a Dirac gap of 25o- The mass term A is modelled by Eq. (|23p . (c) The contour plot shows the density of states; the 
maximum number of states is restricted to five. 

condition E{m, t)——E{1 — m, r) and thus the DOS is symmetric with respect to i? — > —E. For a smah 6o only a 
few energy levels near zero lie within the gap, corresponding to small angular momentum states. The general trend 
is that with increasing Sq the number of discrete levels that falls in the gap increases since gradually larger m values 
give confined states. As an example, for Sq=5 meV the discrete levels correspond to angular momentum numbers 
— 1 2, whereas for Sq—20 meV to —9 10. Still, however, for energies near the middle of the Dirac gap, 

—5 meV<i?< 5 meV, the DOS remains low regardless of Sq, because only small m values contribute to this energy 
range. This suggests that the formation of a gap-induced dot can be probed using similar charge measurements as, 
for example, in GaAs dots. The Fermi energy of the graphene dot has to be adjusted near the middle of the energy 
gap where the DOS is expected to be low and scanning tunneling microscopy could be used to probe the confined 
quantum states. 

V. EFFECT OF A QUANTUM WELL POTENTIAL ON A GAP-INDUCED DOT 

In this section the effect of an electrostatic quantum well potential on a gap-induced dot is examined. It is assumed 
that the potential can be generated by a gate electrode and is not related to the presence of the Dirac gap. We 
consider the most common experimental regime where both the electrostatic potential V and the mass term A are 
constant asymptotically. In this case, if Ag—O then confinement depends on energy, as can be seen from Eq. (jl3[) . and 
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FIG. 5: (a) Energy levels versus the electrostatic potential depth Vb, for m=6, t—+1 and a constant Dirac gap of 2A=50 meV. 
(b) As in (a) but for a spatially modulated Dirac gap with an asymptotic value of 25o=50 meV. (c) As in (b) but for m—2. 
The states of the energies marked by circles are shown in Fig. [6] 

it is achieved when {V — E)^ — A'^<0. If Ag has a power-law dependence the states are always confined independent 
of their energy, but in this work we focus on the case where Ag=0, and therefore there is no magnetic field in the 
graphene sample. 

The interplay of a quantum well potential with the gap-induced dot has a drastic effect on the resulting states. In 
particular, gap-induced dot states couple to the states induced by the potential. As a result the region in which the 
resulting hybridized states are localised is tunable with the strength of the potential. Also, Klein tunneling in the 
electrostatic barrier region occurs. 

It has been experimentally demonstrated that gate electrodes can be suspended above the graphene sheet inducing 
a smooth quantum well or barrier potential that is tunable with the applied gate voltagci^^^ The exact potential 
profile can be determined by solving the Poisson equation within a semiclassical Thomas- Fermi model properly taking 
into account charge screening effects and the specific electrode geometryj^iSI However, a smooth and slowly varying 
potential can be modelled, to a good approximation, by the Gaussian expression^i^ 



The parameters Vg and Iq model the effective depth and width of the quantum well potential, which are controlled 
by the geometry of, and the applied voltage to the gate electrodes. Typical gate voltages of a few volts generate a 
potential well depth of some hundreds of meV, with an effective width of some hundreds of nm. This work is concerned 
only with the regime where the quantum well potential is formed inside the spatial region where A<5oj with A as in 
Eq. ([25)1 . therefore for all the calculations Zo=180 nm and Vq ^ 120 meV. 

The energy level diagram as a function of Vq is shown in Fig.[5{a), for a spatially- independent (constant) Dirac gap 
with 2A=50 meV. Confined states have discrete energy levels that lie within the gap, whereas deconfined states form 
an upper and a lower band of (quasi) continuum of levels, for energies E>A and E<—A, respectively. If Vb is small, 
the angular momentum delocalizes the states; therefore confined states are formed after a critical value of Vb, that 
is ^19 meV in Fig. [Sja). As Vb increases, the number of discrete levels increases due to the stronger confinement. 
Energy levels emerge from the upper continuum into the gap, while the lowest discrete levels merge into the lower 
continuum via anticrossing points. These are formed inside the two bands as shown in Ref. [T^ The quantum states 
undergo a transition from confined to deconfined and vice-versa. Specifically, a confined state transits to a deconfined 
via tunneling into the valence (lower) band, when its energy falls below the gap.— In this regime the deconfined state 
has an oscillatory tail but has a large amplitude near the origin of the dot (see below). 

Figures [5Jb) and (c) show the energy diagram when the Dirac gap is spatially modulated with an asymptotic value 
of 2(5o=50 meV. The appearance of discrete levels for Vo=0 is a consequence of the gap-induced confinement. Such 
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FIG. 6: (a) Quantum states for m=2 and for different electrostatic potential depths Vb; from top to bottom: Vb=0, 30, 50, 
70 meV for the left panels, and Vo=0, 30, 70, 100 meV for the right panels. The Dirac gap is spatially modulated with an 
asymptotic value of 2(5o=50 meV. Left (right) panels show states with energies marked by • (o) in Fig. [Sj^c). (b) As in (a) but 
for m=6 and Vo=0, 40, 85, 100 meV for the left panels, and Vo=0, 40, 85, 120 meV for the right panels. Left (right) panels 
show states with energies marked by • (o) in Fig. [5][b). The vertical axes of the insets range from to 1x10^'^. 

levels are absent in the constant gap system. Furthermore, for Vo=0 there are more discrete levels for m=2 than 
TO=6, as expected based on the above arguments. Anticrossing points, formed now between discrete levels, indicate 
a coupling between potential- and gap-induced states. The gap-induced states of the upper energy ladder couple 
strongly to the potential states and for this reason the corresponding anticrossing points are not well-formed. On the 
other hand, states of the lower energy ladder couple weakly to the electrostatic potential for to=6 [Fig. EJb)], and 
therefore the anticrossing points are well- formed, while this coupling is much stronger for to=2 [Fig. [5jc)]. 

Figure [6] illustrates the dot states as the potential depth Vb increases. Consider first the to=2 states shown in the 
left panels of Fig.[BJa). Confinement for Vb=0 is due to the spatial modulation of the Dirac gap and cannot be realised 
for a constant gap. As Vb increases, the gap-induced state couples to the lowest energy state induced by the potential 
well (Vb=30 meV), and its energy decreases. The state acquires a large amplitude in the potential well region and 
with increasing Vb it couples successively to gap-induced states of the lower energy ladder (Vb=50 meV). In this range 
the resulting hybridised state is spread with significant amplitude over the whole non-asymptotic region defined by 
A < (5o. Further increasing Vb, the energy of the state falls below the gap (Vb=70 meV); the state tunnels in the 
valence band and becomes deconfined with an oscillatory tail. 

Strong coupling between gap- and potential-induced states occurs also for the to=2 states shown in the right panels 
of Fig. [S]^a) ; for instance, for Vb=70 meV and 100 nieV. With increasing Vb, the gap-induced state that corresponds 
to the highest energy in the lower ladder (for Vb=0) couples successively to excited potential-induced states. This 
effect is reflected in the energy level diagram [Fig. [5]Jc)] by the appearance of a series of anticrossing points as Vb 
increases. A similar effect is displayed by gap-induced states of smaller energy, although with increasing Vb the number 
of gap-induced confined states in the lower energy ladder decreases. For instance, as seen in Fig. [5l[c) for Vb=0 there 
are five discrete levels, whereas for Vb=50 meV there are three. The number of states can be determined using a 
similar approximate formalism as that in Ref. [ol The necessary condition for the existence of confined states is that 
there must be a spatial region near r=R in which A:|>0, with k\ as given in Eq. ((T^), whereas asymptotically A:|<0. 

In Fig. Oa) the coupling between potential- and gap-induced states of the lower ladder is strong and therefore the 
resulting hybridised states have a relatively strong amplitude over the whole region where A < (Jq. In contrast, the 
coupling is weak for the m=6 states shown in Fig. [lib). Therefore, for Vb=85 meV, the states peak mainly in the 
region defined either by the potential well (left panel) or the Dirac gap modulation (right panel). These two regions 
have a small overlap when m is large and the effective width of the potential well is smaller than the radius of the 
zero-gap region. 

The m-dependent coupling of the states is related to the width of the classically forbidden region formed between the 
confinement regions defined by the gap modulation and the electrostatic well. The semiclassical approach developed 
in Ref. [H shows that for a large positive to the width of the forbidden region is large. This results in a weak coupling 
between the potential- and gap-induced states of the lower energy ladder. Also, as to increases the gap between the 
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FIG. 7: (Color online) Density of states versus energy E {\E\ < 24.8 meV) and electrostatic potential depth Vo, for r=+l, and 
an asymptotic Dirac gap of 2(5o=50 meV. 

upper and the lower energy ladders increases too, (Fig.U]) and therefore Vo has to be large for the formation of coupled 
states. The coupling effect, though, is insensitive to the details of the gap modulation. Our calculations show that 
both weak and strong couplings can be realised even when the Dirac gap changes discontinuously as in Eq. (jl7p . 

When the potential is nonzero (Vot^O) the induced states can have an oscillatory amplitude inside the barrier region, 
though asymptotically they decay.— This effect is a consequence of the Klein tunneling in the non-asymptotic region 
defined by A < (Jq. In contrast, the quantum states of a Schrodinger dot display only an exponential decay inside an 
electrostatic barrier. Numerical calculations of the DOS, shown in Fig. [3 suggest that it should be experimentally 
possible to resolve graphene states with large oscillatory amplitude in the barrier and the coupling between the states 
described above. In order to probe the states with scanning tunneling microscopy the DOS has to be low. This could 
be achieved for a small Vq, since in this case only energies that correspond to small m states fall in the gap. For the 
same reason the DOS is low when the asymptotic value of the Dirac gap is small. 

VI. DISCUSSION AND CONCLUSIONS 

The existence of a Dirac gap in the energy spectrum of monolayer graphene suppresses the Klein tunneling and 
thus enables charge confinement by an electrostatic quantum well and formation of quantum dots. Here it was shown 
that a graphene dot can be formed as a result of a spatially modulated Dirac gap. It was found, by solving the Dirac 
equation in the continuum limit, that when the gap has a local minimum confined states with discrete energy levels 
can be formed, without applying external electric and/or magnetic fields. This cannot be achieved in a constant Dirac 
gap graphene system. The required gap modulation may be introduced with substrate engineering, local strain, or a 
chemical technique. Unlike quantum dots formed in nanocrystals of graphene, the proposed dot is formed in a large 
sheet of graphene, therefore the physics of the edges is unimportant and does not affect the dot properties. 

The general conditions for confinement in the presence of a spatially modulated Dirac gap, an external electrostatic 
potential, and a magnetic vector potential were analysed. When the gap and/or the vector potential rises asymp- 
totically, while the electrostatic potential is zero (constant), the resulting states are confined regardless of angular 
momentum, valley and eigenenergy. The interplay of the gap-induced dot with an electrostatic potential leads to 
tunable quantum states, i.e., from confined to deconfined and vice- versa. This can only happen provided the potential 
and the mass term, which generates the Dirac gap, have the same spatial forms asymptotically, while the magnetic 
field is zero. In this case the states can be tuned with the strength of the potential: the states are deconfined when 
the electrostatic potential is stronger than the mass term and confined in the opposite regime. 

When the Dirac gap is zero within a disc area and constant outside that area confinement is energy-dependent, 
thus the choice of angular momentum and valley is important. Confined states are localised inside the disc area when 
their energies lie in the gap, otherwise the states are deconfined. The energy spectrum of the confined states consists 
of two ladders (sets) of discrete levels separated by a gap. Numerical calculations of the DOS suggest that states with 
small angular momentum lie in a region of low density and hence they could be probed using standard techniques 
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such as scanning tunneling microscopy. 

It was also shown that states induced by a quantum well potential, which may be generated by a gate electrode, can 
coexist and couple to gap-induced dot states. When the coupling is weak the states retain their character, whereas 
in the opposite regime the states are strongly hybridised. The signature of this coupling in the energy spectrum is 
the appearance of a series of anticrossing points. This coupling property offers a way of tuning the spatial region, in 
which the hybridised states are localized, with the strength of the potential. Moreover, the states can have a large 
oscillatory amplitude in the barrier region exhibiting Klein tunneling, before they eventually decay asymptotically. 
Calculations of the DOS indicate that these quasi-relativistic effects could be probed. 
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